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Abstract 

We consider the acceleration of particles (protons) in a plasma clot by an induced 
field which appears due to the dissolving of the clot. The mechanism of this acceleration 
is based on a model where the clot consists of double flat moving current layers. It 
is assumed that the currents flow in opposite directions, while the clot moves in a 
perpendicular one. It is shown that protons can reach the highest observed energy (of 
the order of 10^^ — lO^^GeV), providing that the size of the clot is large enough. 
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There is a well known problem in high energy cosmic rays physics of how one can explain 
the existence in the universe of protons with the energy in the range of 10^^ — lO^^GeV. 
Unfortunately all existing models do not successfully describe this phenomena. 

In this paper we are going to discuss a mechanism that might explain the existence of 
particles which have such high energies. This mechanism assumes that somewhere in the 
universe there exist very large and sufficiently dense plasma clots. Such objects recently 
were discovered in space, indeed these are the SS433 objects [0] and the jets from the active 
galactic nuclei 0. We also assume that inside the clots there exists an internal current as it 
will be described below. 

Plasma clots usually posses collective electromagnetic field. This field may be produced 
by either a collective electric charge or a current inside the clot. The collision of the clot's 
particles with the medium, as well as the internal recombination in the cooling clot, are 
responsible for dissolving of the moving clot. Therefore, the collective electromagnetic field 
decreases, and the energy and the momentum of the vanishing field can be distributed 
among the rest of the clot particles. This means that the energy of the collective field can be 
transferred to a comparatively small number of particles, and those particles can reach high 
energy. This effect might explain the origin of cosmic ray protons with the highest observed 
energies. Recently, some events were discovered where the energies are larger than 10^^ GeV 



As the first step toward the justification of the above proposed acceleration mechanism 
, we shall show that the energy spectrum of accelerating particles is close to the observed 
one. For this purpose it is natural to assume that the number of particles, say Np, in the 
clot decreases exponentially in time 



where A^o p is the initial number of particles in the clot and a is a slowly varying function of 
time. Then, during the time interval dt the loss in the clot energy, dE, can be presented as 



where e is the sum of the kinetic energy of a particle and a portion of the collective field 
energy per one particle at a given instant of time. A part of energy dE is transferred to the 
rest of the clot particles, such that each of them receives the energy portion de 
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dE = edNp 
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where k is a proper fraction. This gives 
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Then a simple relation holds 
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This spectrum is close to the observed one if the parameter k is of the order of one-half 
to two-thirds. This shows that a /c-portion {k ~ 0.5) of the energy of the clot's collective 
field per one particle can be transferred to the rest of the particles inside the clot. 

To estimate the maximum of the energy the particles (protons) can gain, let us consider a 
simple model as follows. We assume that the collective electromagnetic field is produced by 
a current inside the clot. This current consists of the two double fiat current layers flowing 
in opposite directions, say x-direction. The current inside the single layer is composed by 
both the electrons that have the velocity v^^^ and protons that have the velocity We 
also suppose that both layers as whole system is moving in the perpendicular direction, say 
the z-one, with the velocity Vz relative to the medium. 



Due to the above proposed model the current density in the reference frame K' where 
the current layers are at rest can be expressed in the following form 



f, if 0<z' <L' 
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where L' is the effective width of the layer. 

The corresponding vector-potential A' and magnetic field B' are (we present only the 
none zero components ) 



A' 



B' 



Hoj'z'{L' - z'/2), if 0<z' <L' 
Hojz'lv + z'l2), if - L' < < 
Ho3L'^/2, if z' > L' 
-/xo/L'V2, if z' < -L' 

l^oj'iL' - z'), if 0<z' <L' 
IXoAl' + z'), if -L'<z'<0 
0, if \z'\ > L' 
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Due to the collisions with the medium particles part of the electrons and protons inside 
the current layers live the clot. Because the plasma clot remains quasineutrality both number 
of particles that leave the clot are the same. Therefore current density j'^ decreases. In the 
reference frame K' medium particles have the velocity Vz directed opposite to the z-axis. 
This velocity Vz is the one of the current layer in the laboratory reference frame K where 
the medium is at rest. In the reference frame K' the concentration of the clot electrons n' 
inside the layer can be written in the following form 

n' = n'^expi-aN'vzt') (9) 

where n'^ is the initial concentration of the clot electrons in the layer, A^' is the concentration 
of the particles in the medium, and a is the scattering cross section of an electron inside the 
clot by the medium particles. Here we neglect the thermal motion of the electrons inside the 
clot. 

Let us now pass to the laboratory reference frame K. The component A'^ of the vector 
potential (^ does not transform because it is the transverse component of a 4-vector (as 
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well as the current density j), although the variables it depends upon, do. Then this vector 
potential in the reference frame K has the form 

= ^j{z-v,t)[2LT{z-v,t)H, (10) 
j = en'f^v^expi-aNv^it - ^z)'jI}, (11) 

7. = (1-4)^^ (12) 

where the sign + corresponds to the region < z — Vzt < L, while the sign — corresponds to 
the region —L < z — v^t < 0. The corresponding electric and magnetic fields are (non zero 
components only) 

= ^ya + Bye 

B% = li^3[LT{z-vMll B^c 
= + E^ 

^XC = ^zBy^, E^^ = —By^ 

The fields E^^ and S± are the usual fields produced by the moving current, while the fields 
Ef^ and By^ are the ones which appear due to the dissolving of the clot. The acceleration 
effect we consider here (acceleration of the clot particles) is due to the fields E^^ and B^^. 

We estimate the maximum of the energy the particles can gain due to the described above 
acceleration. Notice, that the potential (p!0| ) does not depend on the transverse coordinate 
X. Therefore there exists an integral of motion 

+ eA^ = pox + eAo^ = const. (15) 

where the index labels the initial condition. As j^O, the potential vanishes as well. 
This allows one to find the maximum of the transverse momentum 
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v,t)['2L±{z-v,t)]x 
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(Px)max = POx + eAo^ (16) 

Let us consider now the longitudinal motion of the particles. The magnetic field -B^ produces 
the Lorentz force 

= ev^B^, (17) 

that acts on the electrons as well as on the protons in the positive z-direction in both regions 
z — Vzt > and z — Vzt < 0. It should be noticed that the field B^^ is responsible for 
the acceleration of the clot as the whole system in 2;-direction. Some evidence of a such 
acceleration were mentioned in 

Notice that the longitudinal force which is acting on the electrons is greater than the 
one acting on the protons. This is because the transverse velocity of the electrons is 
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greater than the one for the protons. Due to the quasineutrahty of the clot, there appears a 
longitudinal electric field Ez that decelerates the electrons and accelerates the protons. 

The longitudinal equations of motion for both electrons (e) and protons (p) are 

§ = -eEf + ev:Bf,, (18) 

As a result, one obtains 

^^^^ = eK + v^.)Bt. (19) 



The longitudinal force in eq.(19) is Vz{vl + v^)/2c^ times less than the force 2eE^ which is 



acting on the e^p pair in the current direction. Therefore the acceleration in the x-direction 
is at least not less than the one in the z-direction. This allows one to estimate the maximum 
of the proton energy 

Cmax ~ (Px)maxC = c{p^O + eA^o) (20) 

If the increase in the energy is much larger than the initial energy, then one finds 

Cmax ~ ecA^ ^ yecjo(2;o - Vzo){2Lo =f {zq - Vzot)hlo (21) 

where jo is the initial current density. 
One can accept an estimation 

{zo - Vzot)i2Lo T (zo - Vzot))-flo ^ L'^ (22) 

This leads to an estimation for the maximum energy 

Cmax ~ yecjo-^^o^T^o (23) 

Thus we can write an estimation for the maximum gained energy 

e^ax ~ 2jo4'7- olO-'(Gev) (24) 
The relatively new evidences and models show that the values 

L'o ^ (0.01 -0.1)pc 

do not appear to be unrealistic. Accepting this estimation we can find the current density 
jo which is necessary for the protons to gain the maximum energy of the order of lO^^Gev. 
This gives 

jo ^ 5(10"" - 10"^2)A/m2 (25) 
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Estimate now the initial magnetic field which corresponds to the current density (pS]). As- 
suming that the relativistic factor 7^ is of the order of one we get from (0) the following 
relation 

Bo ~ /Uoio^o (26) 

This yields 

BoLo ^ ^ (27) 

ec 

For emax ~ lO^^Gev and Lq ^ O.lpc the necessary value of the magnetic field is about 3Gs. 

These estimations do not depend on the model. In order to prove this, we consider two 
different models 

1) Thin long coil. The vector potential which describes the current at the clot particles 
location has the form 

= if jar (28) 
where r is the radius of the coil, and a is its thickness. 

2) Thin round current loop. The vector potential describing the current at the clot 
particles location has the form 

A^ = ^jaHn{8R/a) (29) 

where a is the radius of current string, and R the radius of the loop. 

It can be shown that in these models the result (^) for the maximum energy has a 
similar form if one makes the following substitutions 

L^^ar for the model ( pH]) 
Lg^^a^ ln(8i?/a) for the model (|9D 

Similar effects appear when the collective field inside the clot is produced by an effective 
electric charge in the clot instead of the current. 
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